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In order to investigate the size limit on spatial localized structures in a nonlinear system, wc 
explore the impact of linear nonlocality on their domains of existence and stability. Our system 
of choice is an optical microresonator containing an additional metamaterial layer in the cavity, 
allowing the nonlocal response of the material to become the dominating spatial process. In that 
case, our bifurcation analysis shows that this nonlocality imposes a new limit on the width of 
localized structures going beyond the traditional diffraction limit. 
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Localized structures, belonging to the class of dissipa- 
tive structures found in systems far from equilibrium, arc 
stationary peaks in one or more spatiotemporal dimen- 
sions. They occur in diverse fields of nonlinear science, 
such as chemistry [1, 2], plant ecology [3], gas discharge 
systems [4], fluids [5] and optics [6-12]. The presence 
of localized structures has been demonstrated in systems 
exhibiting a nonlocal response, such as in models for pop- 
ulation dynamics [13], neural networks [14] and nonlocal 
optical materials [15, 16]. Although the formation and 
stability of these localized structures must be influenced 
by these nonlocal effects, the main effects of nonlocal- 
ity are typically obscured by the presence of a stronger 
diffusion and/or diffraction process. Thanks to recent 
advances in the fabrication of mctamatcrials, it is now 
possible to conceive a nonlinear system under conditions 
such that the dynamical properties of localized structures 
are dominated by the nonlocal response. In this work, it 
is our aim to analyze the influence of this nonlocality on 
the nonlinear dynamical properties of two-dimensional 
(2D) spatial structures. 

Localized structures are relatively well understood in 
one transverse dimension [17], whereas an analytical 
analysis in two transverse dimensions is still largely un- 
explored and most of the results are obtained by numer- 
ical simulations. In this work, we will focus on the par- 
ticular system of controllable 2D localized structures in 
optical cavities, also referred to as cavity solitons (ess). 
Their formation can be attributed to the balance between 
nonlinearities due to light-matter interaction, transport 
processes (diffusion and/or diffraction), and dissipation 
[7, 8, 18, 19]. These bright spots have been proposed 
for information encoding and processing [8, 18, 20]. De- 
creasing the size of CSs would be advantageous for these 
applications, while also being of fundamental interest. In 
optical devices, the spatial dimension of CSs is generally 
restricted by diffraction, which imposes the diameter to 



be of the order of \/V, where V = IT j (nk) is the diffrac- 
tion coefficient, with k the wavenumber of the beam, and 
I and T the length and the finesse of the resonator. For 
example, in Ref. [8], CSs at 0.5|J.m have a transverse size 
that reaches the diffraction limit of 10 (xm. 

The use of transverse index modulation has been pro- 
posed to overcome the diffraction limit of mid-band CSs 
in a certain class of resonators [21] . These localized struc- 
tures have a large intrinsic transverse velocity that pre- 
vents its use for many practical applications. Recently, 
a different strategy was motivated by the progression of 
the fabrication of left-handed materials (lhm) towards 
optical frequencies [22, 23], and the proposition of non- 
linear lhm [24, 25]. In Rcfs. [26, 27], Kockaert et al. 
study the possibility of altering the strength of diffrac- 
tion by inserting a layer of lhm, in addition to a layer of 
right-handed material (rhm), in the cavity of an optical 
microresonator, and they show how to reduce the diffrac- 
tion coefficient to arbitrarily small values. As the soliton 
width scales with the square root of the diffraction coeffi- 
cient, this method potentially allows for sub-diffraction- 
limited CSs. Unlike with the use of photonic crystals, 
this technique works in principle for all types of CSs and 
microresonators. Although the diffraction limit can be 
encompassed in this way, one can reasonably expect that 
the sub-wavelength structure of lhms will impose a new 
size limit on CSs. Here, we will show that the inherent 
nonlocality of lhms will significantly change the proper- 
ties of CSs when diffraction is tuned down, and a new 
size limit of CSs — now imposed by nonlocality — will 
be revealed. 

We consider a microresonator driven by a coherent op- 
tical beam. In each roundtrip, the light passes through 
two adjacent nonlinear Kerr media: a RHM and a LHM. 
It has been shown in Ref. [26] that the evolution of the 
electric field in this microresonator is governed by the 
well-known Lugiato-Lefever (LL) equation [28], with the 



diffraction coefficient V given by 

2n 2 Ur K\J [ ' 

1l ! an d II are the indices of refraction and the 
lengths of the RHM and the lhm, respectively. By chang- 
ing l R and 2? can be engineered to ever smaller, pos- 
itive diffraction coefficients. From the LL equation, one 
can estimate the diameter of localized solutions to be 

with A the cavity detuning. Therefore, CSs will become 
infinitely small when T> tends to zero. But in that case, 
higher order effects will start to dominate the spatial dy- 
namics. Indeed, from the derivation of the LL equation, 
one can show that the inherent nonlocality of the nanos- 
tructured metamaterial comes into play when diffraction 
becomes negligible. The nonlocality in these materials is 
still largely unexplored. When addressing the cavity with 
an optical beam, the induced fields in the sub-wavelength 
resonators will also couple to neighboring resonators, pro- 
viding a nonlocal response. We assume this nonlocal re- 
sponse to be weak, because the coupling to the nearest 
neighbor resonators is important. 

Under the same approximations under which the LL 
equation is valid, i.e., slowly varying envelope approxi- 
mation, weak nonlinearity and a nearly resonant cavity, 
we find that the nonlocality comes in as follows: 

8 A 

— = -(1 + iA)A + A in + i \A\ 2 A + WVlA 

0(t ± - u)A(u)du. (3) 

The kernel function 8, describing the nonlocal response 
of the linear lhm, effectively couples the electric field 
at different positions. The details of the derivation of 
Eq. (3) will be published elsewhere. When the nonlocal- 
ity is weak, the last term of Eq. (3) can be expanded in a 
series of spatial derivatives of A, and, taking into account 
the rotational invariance of the system, we find 

J 1 6(r ± - u)A(u)du ~ 9 A + 9 1 \/ 2 L A + 9 2 \/\A. (4) 

The first two terms in Eq. (4) only change the value of 
the detuning and diffraction coefficients, respectively. In 
what follows, we absorb their contribution into the pa- 
rameters A and keeping the same notation. We 
finally arrive at 

^ = -(l+iA)A+A in +i\A\ 2 A+W^VlA+W^ViA. 
at 

(5) 

Eq. (5) is similar to the LL equation, and has the 
same homogeneous steady state (hss) solutions ^4 S . In 
this work, we want to study the localized structures 
arising from the modulational instability of the homo- 
geneous solution. Therefore, we restrict ourselves to the 



2 




FIG. 1. Influence of the nonlocality to diffraction ratio 77 
on the marginal stability curves as given by Eq. (6). The 
wavevectors inside the curves destabilize the homogeneous 
steady state solution of Eq. (5). A = 1.23. (a) Positive 77. 
(b) Negative 77. 

monostable regime (A < \/3). We have performed a sta- 
bility analysis by linearizing Eq. (5) around the steady 
state solution, and seeking for the deviation in the form 
A = A s + 5A exp(ik-r + At), with k = (k x ,k y ) and 
r = (x,y). The marginal stability is given by 

V^k 2 m - V^k 4 m = 2\A S \ 2 — A ± Vl^l 4 - I- (6) 

The modulational instability (MI) depends strongly on 
the parameter 77, defined as the nonlocality to diffraction 
ratio rj = 2?( 2 ) /V^ . In Fig. 1, the marginal stability 
curves are shown for several values of 77. 

For 77 > [see Fig. 1(a)], these curves have the form 
of a cardioid, which contains the unstable wavevectors. 
A first MI arises with two critical wavenumbers at an 
intensity of \A S \ 2 = 1. This leads to the formation of 
complex patterns with two wavevector components. Sta- 
bility of the HSS is recovered at higher background inten- 
sities. As r\ increases, the cardioid form evolves into an 
elliptical shape with equal critical wavenumbers at both 
sides: T>^k^ = 1/47/. From this wavenumber, we can 
estimate the typical width of localized structures to be 
of the order of 

A+ = 27rv / 4£>( 2 )?7. (7) 

In the limit of small diffraction, Eq. (7) predicts that the 
width of CSs increases when T>^ decreases, while Eq. (2) 
indicates a decreasing width. Therefore, there must exist 
a value of the diffraction strength for which the CS 
is of minimal width. We will use a numerical method 
to determine this optimum. However, we want to point 
out that this minimal size CS can be unstable. Our cal- 
culations will also provide insight in this matter. Also 
note that, when 77 > 1/4 (\/3 — A), the unstable region 
contained in the marginal stability curve disappears, ef- 
fectively stabilizing the HSS. 
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FIG. 2. Maximal intensity of the cavity soliton vs. the back- 
ground intensity for different values of 77. Solid lines indicate 
stable structures, whereas dashed lines correspond to unstable 
solutions. A = 1.23. (a) Positive r\. (b) Negative 77. 

For r] < [see Fig. 1(b)], the modulational instability 
occurs again for intensities |j4 s | 2 > 1, but here the HSS is 
not stabilized for high input fields. At the critical point, 
we can derive the following estimate for the CS width: 

A = , 2 -^^ . (8) 
yJl + y/l-4(2-A)r, 

Consequently, when 77 < 0, the CS width will still 
decrease with the diffraction strength but due 

to the nonlocality the width will saturate at Ar lm — 
2tt{/-T>( 2 )/(2 - A). Again, the possibility that the CSs 
become unstable before reaching this limit exists, empha- 
sizing the need to check the stability of these structures. 

Our numerical analysis relies on a calculation of the 
stationary localized solutions of Eq. (5) using a Newton- 
Raphson method [29-31]. The radial form of this equa- 
tion is discretized, from which a set of coupled nonlinear 
equations is obtained. Since the equation is linear in 
the spatial derivatives, the diffraction term can be com- 
puted in the spatial Fourier space. Zero derivatives at 
the boundaries are imposed. This approach is extremely 
accurate and generates the Jacobian. 

In Fig. 2, CS bifurcation diagrams are shown for differ- 
ent values of 77. In the absence of nonlocality (77 — 0), as 
previously studied in Refs. [29, 30], a CS branch emerges 
subcritically from the HSS at |^4 S | 2 = 1. The negative 
slope part of this branch corresponds to unstable CSs; 
the positive slope higher branch is stable for low val- 
ues of I A s 1 2 . At higher background intensities the in- 
tensity peaks become higher and narrower, attributed to 
the more prominent self-focusing effect, and at a certain 
point they become again unstable with respect to az- 
imuthal perturbations. For 77 > [Fig. 2(a)], we observe 
that the branches extend to smaller |^4 S | 2 with increas- 
ing 77, while the domain of stability is enlarged. This 
trend is reversed from a certain 77 on, and the stability 



FIG. 3. Renormalized full width at half maximum £ vs. 
diffraction strength (M^ 1 / 2 ). The width £ is shown of the 
higher branch CS at the change of stability or of the CS at the 
saddle-node bifurcation when no stable CSs exist Stable CSs 
are indicated by solid lines, whereas the unstable CSs are in 
dotted lines. A = 1.23. (a) Positive 77. (b) Negative rj. 

range decreases. Also note that due to nonlocality, which 
physically tends to spread out the inensity to neighbor- 
ing points, the peak intensity of CSs drops. For rj < 
[Fig. 2(b)], the branches move monotonically to higher 
background intensity and to lower peak intensity with 
stronger nonlocality. Again note that the stability range 
is reduced considerably. 

To find the minimal CS size as discussed above, we have 
investigated the scaling of the renormalized full width at 
half maximum (£ = dFWHM / \/|X>( 2 )|) of the CSs with rj 
in Fig. 3. Note that this renormalization enables us to 
obtain a general result only depending on the detuning 
A. For each value of 77, we indicate the width £ of the 
upper branch CSs at the change of stability, when stable 
CSs exist, while plotting the £ at the saddle-node bifur- 
cation when the entire upper branch is unstable. For 
77 > [Fig. 3(a)], one can distinguish two regions of sta- 
ble CSs, one in which the width decreases with decreasing 
diffraction strength, and the other with the inverse effect. 
One can identify the minimal width £ min and its corre- 
sponding optimal 7/- op t in the leftmost region. For rj < 
[Fig. 3(b)], the width decreases monotonically with de- 
creasing diffraction strength, until stability of the CSs is 
lost. At this point, the minimal width £ mirl is obtained. 
Our numerical calculations thus qualitatively confirm the 
prediction by the linear stability analysis given above. 

We have repeated this procedure for several values of 
the detuning. This gives us the optimal parameter set 
that produces minimal size CSs given a certain nonlo- 
cality X>( 2 '. Fig. 4 summarizes the main results. For 
77 > [Fig. 4(a)], one can see that the minimal width 
£min varies only by a few percent in this range of the 
detuning. The smallest £ m ; n can be obtained for higher 
values of the detuning A, and this with smaller values 
of rj opt . For the corresponding optimal values of rj, 
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FIG. 4. Minimal CS width £ m i n (solid lines, left axis) and opti- 
mal rjopt (dashed lines, right axis) vs. detuning A. (a) Positive 
t/. (b) Negative 77. 

remains smaller than V^ 2 , so the formation of the min- 
imal width CSs is still dominated by diffraction effects. 
For r/ < [Fig. 4(b)], the minimal width again decreases 
with increasing values of the detuning, but changes more 
strongly than for positive r/. Note that this time absolute 
values of rj opt need to be larger to obtain the smallest CSs. 



Hence, for r\ < 0, it is the nonlocality that has the largest 
influence on the formation of minimal size CSs. 

In conclusion, we have studied the impact of non- 
locality on the nonlinear dynamical properties of two- 
dimensional spatially localized structures. Our study 
confirms the possibility to reduce the size of cavity soli- 
tons beyond the diffraction limit by using left-handed 
materials. However, we have shown that the nonlocal in- 
teraction between the optical field and the material not 
only hinders this size reduction, but also alters the stabil- 
ity of cavity solitons in a manifest way, imposing a new 
limit on their size. 
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